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Vortex critical dynamics of the two dimensional XY spin glass is studied by Monte Carlo methods in 
the Coulomb-gas representation. A scaling analysis of the nonlinear response is used to calculate the 
correlation length exponent v of the zero-temperature glass transition. The estimate, v — 1.3(2), 
is in agreement with a recent estimate in the phase representation using the same analysis and 
indicates that the relevant length scale for vortex motion is set by the spin-glass correlation length 
and that spin and chiralities may order with different correlation length exponents. 
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It is well known that vector spin glasses, such as the XY 
spin glass, have a chirality order parameter with Ising-like 
symmetry in addition to the continuous degeneracy as- 
sociated with global spin rotationtl. Chirality arises from 
quenched in vortices due to frustration effects in each el- 
ementary cell of the lattice which contains and odd num- 
ber of antiferromagnetic bonds. The interplay between 
spin and chiral variables has always received considerable 
attention because of the possibility of separate spin-glass 
and chiral-glass ordering djLie_to freezing of spins and chi- 
ral variables, respectivelyQ'Erffl'El. Possible separation of 
spin and chiral variables also arise in the frustrated XY 
model with weak disorderLla. While in three dimensions 
the existence of a fiiM^A tepiperature transition is under 
current investigationEi'Eja, in two dimensions there is a 
consensus that the transition only occurs at zero temper- 
ature. Associated with the zero temperature transition 
there is a correlation length which increases with decreas- 
ing temperatures as ^ oc T^". However, the possibility 
of different spin and chiral glass short range correlation 
lengths S,s and S,c with different critical exponents Vs and 
i/c has not been resolved satisfactorily. 

First evidence that spin and chiral glass correlation 
length exponents are different in two |dimensions were re- 
ported by Kawamura and Tanemuraa from domain wall 
calculations. Varioua-icstimates of these exponents give 
approximate valuesBuO'El of j^s ^ 1 and i^c '^ 2 but the 
errorbars are usually quite large and a single exponent 
scenario may not be ruled out which would be consistent 
with an analjytical work for a particular type oixlisorder 
distributions and more recent numerical workll3 on do- 
main wall scaling behavior at zero temperature. These 
calculations are usually performed in a representation of 
the XY spin glass model in terms of the orientational 
angle of the two-component XY spins. In this represen- 
tation, spin glass order can be directly identified as a 
long range order in appropriate phase correlations while 
the chiral variables are built from nearest neighbor phase 
correlations. In numerical simulations, the dynamics of 
the chiral variables are then determined by the phase 



variables and equilibration problems may prevent an ad- 
equate study of the vortex correlations in the system. 

An alternative approach which allows to study the vor- 
tex dynamics directiy can be obtained from the CouJomb- 
gas representationllJ. Recently, Bokil and Youngp used 
Monte Carlo simulations in the vortex representation to 
obtain an estimate of the chiral glass correlation length 
exponent and found i^c — 1-8 ± 0.3. This agrees with 
previous estimates within the errorbar. It also supports 
the scenario in which spin and chiral glass variables order 
with different critical exponents if one accepts earlier es- 
timates of the spin glass exponent ^'c ~ 1, obtained in the 
phase representation. However, it should also be of inter- 
est a determination of the spin glass correlation length 
exponent from simulations in the vortex representation 
as it is not completely clear how such length scale shows 
up in the dynamical behavior of vortices. In particular, 
since the XY spin glass model is currently been used as 
a model for gcaiuilar high- Tc superconductors containing 



TT junction: 



which leads to quenched in vortices 



even in the absence of external magnetic field, a natural 
question arises as to which correlation length, ^s or ^c, 
is actually probed by transport measurements. In the 
measurements, the response of the vortices to an applied 
force can be observed as the voltage response to an ap- 
plied drivin g . cu rrent which acts as a Lorentz force on 
the vorticcsOlla. The vortex response, or resistive be- 
havior, is therefore determined by vortex mobility and 
the current-voltage scaling is expected-.to be controlled 
by the relevant divergent length scaleE3 which could be 
either ^s or ^c- 

The question of the relevant correlation length for 
vortex response is also of interest for the three dimen- 
sional XY_spin glass. Recent simulations of the vortex 
dynamicsty in three dimensions showed evidence of a re- 
sistive phase transition at finite temperatures which was 
attributed to glass ordering of chiralities while the spins 
remain disordered. This would be consistent with the 
scenario of a finite temperature chiral glass transition in 
the absence of spin glass transition which has been pro- 



posed previouslycl from calculations in the phase repre- 
sentation of the XY spin glass model. This interpretation 
however is only justified if the relevant length scale for 
vortex dynamics is determined by the chiral glass correla- 
tion length ^c- On the other hand, since it is well known 
that vortex motion leads to phase incoherence, one ex- 
pects that vortex dynamics should probe the spin glass 
correlation length ^^ and therefore the resistive transition 
should correspond instead to a spin glass transition at fi- 
nite temperatures. In fact, this is supported by a more 
recent domain-wall calculations suggesting a spin glass 
transition at finite temperatures in three dimensionstil. 
The present study of the vortex response in two dimen- 
sions may help to clarify this point as it is well known 
that in this case both spin and chiral glass ordering only 
occurs at zero temperature and so the scaling analysis 
involve less unknown parameters. 

In the absence of a precise agreement among the var- 
ious studies of the XY spin glass model, as mentioned 
above, and in view of its relevance for vortex dynam- 
ics, the additional numerical results presented below may 
help to settle some issues. 

In this work, we study the vortex critical dynamics 
of the two dimensional XY spin glass by M«nte Carlo 
methods in the Coulomb-gas representationll3. A scal- 
ing analysis of the nonlinear response is used to calculate 
the correlation length exponent v of the zero-temperature 
glass transition. The estimate, v = 1.3(2), is in agrefe 
ment with a recent estimate in the phase representationEZI 
using the same analysis and indicates that the relevant 
length scale for vortex motion is set by the spin-glass cor- 
relation length ^s and that spin and chiralities may order 
with different correlation length exponents. 

We consider the XY spin glass on a square lattice de- 
fined by the Hamiltonian 
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(1) 



where 9i is the phase of a two-component classical spin 
of unit length, Si — (cos 6*^, sin 0^), J > is a coupling 
constant and Aij has a binary distribution, or tt , with 
equal probability, corresponding to a coupling constant 



'Ji', 



-J or J, respectively, between Si and Sj spins. 



The sum is over all nearest-neighbor pairs. This Hamilto- 
nian also describes an array of Josephson junctions where 
there is a phase shift of tt across a fraction of the junc.tJQns 
as in models of d-wave ceramic superconductoraJ'tirEj. 

To study the vortex dynamics it is convenient to 
rewrite the above Hamiltonian in the Coulomb-gas rep- 
resentation 



Hcg = 2Tr^ J ^{Ur - ,fr)G'„,{nr> - fr 



(2) 



which can be obtained following a standard procedureta 
in which Eq. (1) is replaced by a periodic Gaussian model 
separating spin- wave and vortex variables. The Coulomb 
interaction is given by G^^, = G{r — r') — G(0), where G 
is the lattice Green's function. 



^w-T^E 



exp(z/c • r) 
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(3) 



and L is the system size. G'{r) diverges as 27rlog|r| at 
large separations r. The vortices are represented by in- 
teger charges Ur at the sites r of the dual lattice and the 
frustration effects of Aij by quenched random charges fr 
given by the directed sum of Aij around the plaquette, 
fr — J2 Aij/2TT. The charges are constrained by the neu- 
trality condition J2ri''^r — fr) = 0. For the XY spin glass, 
the charges fr have a correlated random distribution of 
integer and half integer values. Other random distribu- 
tions can represent different models. Ajiiiiform distri- 
bution describes the gauge glass modelHrEj where Aij in 
Eq. (1) is a continuous variable in the range [0, 27r] and an 
uncorrelated continuous distribution can describe arrays 
of superconducting grains with random fluxCI or arrays^ 
mesoscopic metallic grains with random offset chargealZI. 

We study the nonequlibrium response of the vortices 
in the XY spin glass by Monte Carlo simulations of the 
Coulomb gas under an applied electric fielcll3. An electric 
field E represents an applied force acting on the vortices 
and gives an additional contribution to the energy in Eq. 
(2) as ^^ E rir Xr for E in the_jE direction. A finite 
E sets an additional length scaldl^ in the problem since 
thermal fluctuations alone, of typical energy kT , leads to 
a characteristic length I ^ kT/E over which single charge 
motion is possible. Thus, increasing E will probe smaller 
length scales. Crossover effects are then expected when I 
is of the order of the relevant correlation length for inde- 
pendent charge motion. As vortex motion leads to phase 
incoherence we thus expect that the scaling behavior of 
the nonlinear response will probe the spin glass correla- 
tion length ^a of the original model and allow an estimate 
of the thermal critical exponent Vs ■ This dynamical ap- 
proach complements previous equilibrium calculations in 
the vortex representation of the XY spin glass jwhere only 
the chiral glass correlation length was studiedj. 

In the dynamical simulations, the Monte Carlo time 
is identified as the real time t and we take the unit of 
time dt — 1 corresponding to a complete Monte Carlo 
pass through the lattice. A Monte Carlo step consists 
of adding a dipole of unit charges and unit length to a 
nearest neighbor charge pair (n^, Uj), using the Metropo- 
lis algorithm. Choosing a nearest-neighbor pair i,j at 
random, the step consists in changing rii — > n^ — 1 



and 



1, corresponding to the motion of a 



charge by a unit length from i to j. If the change in 
energy is AC/, the move is accepted with probability 
min{l, exp(— AC//A:T)}. The external electric field E bi- 
ases the added dipole, leading to a current / as the net 
flow of charges in the direction of the electric field if the 
charges are mobile. The current / is calculated as 



m 



yE^Q'W 



(4) 



after each Monte Carlo pass through the lattice, where 
L is lattice size and AQi{t) = 1 if a charge at site i 



moves one lattice spacing in the direction of the field E 
at time t, AQi(i) = —1 if it moves in the opposite di- 
rection and AQi{t) = otherwise. Periodic boundary 
conditions are used. Most calculations were performed 
for L — 32 and compared to a smaller system of i = 16 
but size dependence was not significant in the tempera- 
ture range studied. The current density J is defined as 
J = I / L. The linear response is given by the linear con- 
ductance Gl = lim^j^o J /E which can be obtained from 
the fluctuation-dissipation relation as 



Gr = 



1 



2kT 



dt < /(0)/(i) > 



(5) 



without imposing the external field E. In the calcula- 
tions, the integral is replaced by a sum of successive 
Monte Carlo sweeps through the lattice with the unit 
of time dt = 1. We use typically 4 x 10"* Monte Carlo 
steps to compute averages and 20 different realizations of 
disorder. 
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FIG. 1. (a) Nonlinear conductance J/E as a function of 
temperature T. (b) Arrhenius plot for the temperature de- 
pendence of the linear conductance Gl- 

To analyze the numerical results we need a scaling the- 
ory of the nonlinear response near a second-order phase 
transition. A detailed scaling theory has been described 
in the context of the current-voltage characteristics of 
vortex-glass modelsEj but it can be directly applied to 
the present case. Since the glass transition occurs at 
T — with a power-law divergent correlation length 
^ oc T^'^ and the external field introduces an additional 



length scale I oc kT/E, the dimensionless ratio E/JGl 
can be cast into a simple scaling forn£3 in terms of the 
dimensionless argument ^/l, 



J/EGl ^ F{E/T^+'') 



(6) 



where F is a scaling function with F{Q) = 1. This scal- 
ing form indicates that a crossover from linear behavior, 
when F{x) ~ 1, to nonlinear behavior, when F{x) >> 1, 
is expected to occur when a; ~ 1 which leads to a charac- 
teristic field Ec oc T^'^^ at which nonlinear behavior sets 
in. 

The nonlinear response J/E and an Arrhenius plot of 
the linear conductance Gl are shown in Fig. 1. The data 
shows the expected behavior for a T = transition. The 
ratio J/E in Fig. 1(a) tends to a finite value for small E, 
corresponding to the linear conductance Gl in Fig. 1(b) 
with an activated behavior. This activated behavior is 
consistent with a zero temperature transition and finite 
correlation length at nonzero temperatures which leads 
to a finite energy barrier U for vortex motion. In general, 
an energy barrier exponentEEl ip can also be defined from 
[/ ~ ^'^ for a temperature dependent energy barrier. The 
pure Arrhenius activated behavior in Fig. 1(b) is consis- 
tent with an exponent ip ^ 0. As can be seen from Fig. 
1(a), there is a smooth crossover from linear behavior, 
when J/E is roughly a constant, to nonlinear behavior 
for increasing E at each temperature which appears at 
smaller E for decreasing temperatures in agreement with 
the expected crossover behavior at a characteristic field 

Ec oc ri+'^ . 

We now verify the scaling hypothesis of Eq. (6) and 
obtain a numerical estimate of the thermal correlation 
length exponent zy using two different methods. Fig. 2(a) 
shows the temperature dependence of Ec defined as the 
value of E where E/JGl starts to deviate from a fixed 
value of 2. From the expected power-law behavior for 
the crossover field Ec oc T^+'^ we obtain a direct esti- 
mate of 1/ = 1.4(2) in a log-log plot. From a scaling plot 
of the nonlinear response according to Eq. (6), v can also 
be obtained by adjusting its value so that the best data 
collapse is obtained as shown in Fig. 2(b). The data 
collapse supports the scaling behavior and provides an 
independent estimate of v = 1.3. From the two indepen- 
dent estimates we finally obtain i/ = 1.35 ± 0.2. 

Our estimate of i^ = 1.35 ± 0.2 from the scaling anal- 
ysis of the vortex response is consistent with previous 
estimates of the spin glass correlation length exponent 
h's oVjifHiod in the phase representation of the XY spin 
glassaU'B. These calculations give numerical estimates 
with comparable uncertainties ranging from z^s = 1 to 
1.2. It also agrees with a recent calculation in the phase 
representation o£ Eq. {1), i/ — 1.1 ±0.2, using the same 
scaling analysisO. This suggests that the relevant length 
scale for vortex dynamics is set by the spin glass cor- 
relation length which determines short range phase co- 
herence. Since the chiral glass correlation length expo- 
nent 
rangi 
in which the phase and chiral variables in the XY spin 



n estimated to be significantly larger, in the 
Vc — 1.8 to 2.6, it also supports the scenario 



glass are decoupled on large length scales and order with 
different correlation length exponents. However, as the 
errorbars of these, pitimates are significant large, a single 
critical exponentu'Ej may not be completely ruled out on 
pure numerical grounds and further work will be neces- 
sary to completely settle this issue. 
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the effective Hamiltonian of frustrated XY models with 
weak disorderll For a stable decoupled fixed point, the 
coupling term should be an irrelevant perturbation, cor- 
responding to an eigenvalue A = 2 — a: < evaluated 
at the unperturbed fixed point, where 2x is the corrjei 
lation function exponent. Using the scaling relationEj 
X = 2 — l/v for the energy density correlations, and 
the proposed numerical values for the correlation length 
exponents Vs — 1.3 and Vc — 2, we find indeed that 
A = 2 — a;s — Xc < as required for a stable decoupled 
fixed point. If the transition at T = corresponds to a 
decoupled fixed point then phase and chiral variables can 
order with different correlation length exponents. These 
arguments, by no means, show that a decoupled transi- 
tion is actually realized in the XY spin glass but it makes 
plausible the assumption of distinct divergent correlation 
lengths at the same transition temperature used in our 
analysis of the numerical data. 

Finally, our calculation for the two-dimensional XY 
spin glass, which indicates that vortex dynamics probe 
mainly the spin glass correlation length rather than the 
chiral glass correlation length, also suggests that the fi- 
nite temperature resistive transition observed recently by 
Wengel and Young)i3 in numerical simulations in the vor- 
tex representation of the tree-dimensional XY spin glass 
model should be attributed to spin glass ordering. |-J"his 
is in fact consistent with more recent calculation^ in- 
dicating that the lower critical dimension for spin-glass 
ordering may be just above 3. 

This work was supported by Fundagao de Amparo 
a Pesquisa do Estado de Sao Paulo, FAPESP (Proc. 
99/02532-0). 



FIG. 2. (a) Crossover field Ec as a function of temperature, 
(b) Scaling plot J/EGl x E jT^^" for v = 1.3. 



It should be noted that the decoupled scenario for spin 
and chiral variables near the same transition tempera- 
ture, suggested by our numerical results, does not contrau 
diet general arguments of renormalization-group theorjt3 
which allows for the possibility of nontrivially decoupled 
fixed points. In fact, since the model has continuous and 
Ising-like symmetries, one would expect that the effec- 
tive Hamiltonian describing the critical behavior could be 
written in terms of a disordered ferromagnetic XY-spin 
model, with a zero-temperature transition, coupled to an 
Ising spin glass model, representing phase coherence and 
chiral degrees of freedom, respectively. Ferromagnetic 
XY spin models with a zero temperature transition do 
exist as, for example, in diluted XY models at percola- 
tion threshold23'Ell. Although the exact form of the effec- 
tive XY and Ising Hamiltonians and the coupling term 
are not known, the continuous and discrete symmetries 
of the model are consistent with an energy density cou- 
pling of the form J2r^'i(^)^c{r), where Eg and Ec are 
the local energy densities for the XY spins and chirality, 
respectively. Such a coupling term is known to occur in 
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